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INTRODUCTION 

Aerodynamic  ground  test  facilities  for  simulating 
re-entry  flows  involve  the  use  of  very  high  operating 
pressures  at  high  temperatures  (Ref.l).  Under  these  conditions 
the  equation  of  state  departs  markedly  from  the  perfect 
gaB  law  due  to  compres s ihility  effects  and  excitation  of 
internal  energy  modes.  The  associated  thermodynamic  properties 
are  hence  also  strongly  influenced.  Optimisation  studies 
of  ground  test  facilities  require  that  such  gas  imperfections 
are  known  with  a reasonable  degree  of  accuracy.  This  report 
describes  an  experimental  program  to  make  accurate  measur- 
ements of  pressure,  volume  and  temperature  of  nitrogen  and 
air  at  pressures  up  to  3000  kg/cm2  and  temperatures 

up  to  2500°K,a  range  of  conditions  in  which  the  compressibility 
effects  are  most  dominant. 

The  behaviour  of  dense  gases  may,  in  principle, 
be  predicted  using  quantum-statistical  mechanics.  However, 
the  mathematical  difficulties  involved  when  molecular  inter- 
actions must  be  taken  into  account  are  so  formidable  that 
solutions  have  only  been  found  using  highly  simplified  mole- 
cular models  (Ref. 2).  Many  attempts  have  been  made  to  obtain 
semi-empirical  equations  of  state  which  are  valid  for  any 
dense  gas  by  the  application  of  Van  der  Waals ' principle 
of  corresponding  states  (Refs  3 and  4).  These  depend  on 
having  available  knowledge  of  the  equations  of  state  over 
the  range  required.  Once  the  equation  of  state  is  defined, 
all  thermodynamic  properties  can  be  obtained  from  existing 
relations  (Ref. 4). 

Up  to  present  experiments,  providing  the  knowledge 
on  which  the  equations  of  state  are  based,  have  been  confined 
to  the  studies  of  gases  at  high  pressures  and  ordinary 


5 


AEDC-TH-74-47 


temperatures,  where  measurements  may  be  made  under  steady 
state  conditions.  A large  amount  of  data  is  available  in 
this  regime.  The  contrary  is  true  at  high  temperatures 
because  of  the  problems  encountered  in  maintaining  the 
structural  integrity  of  the  containing  vessel.  These  problems 
can  be  overcome  by  rapidly  heating  and  compressing  a gas 
sample  by  a piston,  and  taking  the  desired  measurements 
under  transient  conditions.  Work  of  this  kind  has  been 
carried  out  in  the  USSR,  but  aside  from  a brief  monograph 
on  the  experimental  equipment  and  techniques  employed 
(Ref. 5),  no  results  seem  to  have  been  subsequently  published 
in  the  open  literature. 

The  driver  section  of  the  VKI  Piston  Driven  Shock 
Tube  has  been  modified  to  generate  samples  of  dense  high 
temperature  gases  and  instrumentation  has  been  developed  to 
measure  pressure,  temperature  and  density  (Ref. 6).  A sodium 
spectral  line  reversal  technique  is  used  to  measure  the 
temperature  of  the  gas.  The  density  of  the  gas  is  determined 
indirectly,  by  measuring  the  piston's  position  during  the 
compression  process.  This  measurement  involves  the  use  of 
an  eddy  current  displacement  transducer.  A quartz-crystal 
piezo  electric  transducer  is  used  to  measure  the  pressure 
of  the  gas.  More  details  of  the  instrumentation  are  given 
within.  A full  description  of  the  VKI  Piston  Driven  Shock 
Tube  and  its  control  equipment  before  conversion  to  a 
compression  tube  is  given  in  Ref. 7- 

The  report  describes  experiments  carried  out  to 
generate  thermodynamic  state  data  of  nitrogen  and  air  and 
to  assess  the  experimental  accuracy.  The  results  are  compared 
with  published  thermodynamic  tables. 
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2.  EXPERIMENTAL  TECHNIQUE 

2 . 1 The  test  gas  compression  system 

The  system  for  compressing  the  gas  is  illustrated 

in  Fig.1.  A piston  constructed  of  hardened  steel,  aluminium 

and  nylon,  veighing  12.6  kg  and  of  length  250  mm  is  free 

to  move  in  a barrel  of  12.6  mm  diameter  and  1.9dm  in  length. 

The  barrel  is  attached  to  a reservoir,  of  diameter  270  mm 

and  length  2.3m  at  one  end  and  is  closed  at  the  other  end. 

Initially  an  aluminium  diaphragm,  typioally  of  thickness 

0.4  mm,  retains  the  piston  at  the  reservoir  end  of  the  barrel. 

The  gas  under  test  is  introduced  into  the  barrel  after 

. * 

evacuating  and  purging  the  barrel  with  the  use  of  a vaccum  pump  , 
and  the  reservoir  1b  charged  with  air  up  to  its  operating 
pressure  Pq  at  temperature  Tq.  To  release  the  piston  at  a 
selected  time,  the  reservoir  gas  is  bled  in  such  a vay  as  to 
act  on  a larger  area  of  the  piston  rear  surface,  the  extra 
force  generated  being  adequate  to  shear  the  retaining 
diaphragm.  Subsequently  the  piston  is  driven  down  the  barrel 
compressing  the  test  gas  from  an  initial  pressure  P4-  and 
Tl*i  to  a maximum  pressure  P4f  and  temperature  T4f.  Both  Tq 
and  T4 are  equal  to  the  laboratory  temperature.  The  piston 
is  instantaneously  at  rest  when  the  peak  pressure  is  achieved. 
After  the  first  compression  stroke,  the  piston  reverses  its 
motion  and  returns  to  the  reservoir  end  of  the  barrel.  It 
continues  to  oscillate  in  the  barrel  until  brought  to  rest 
by  friction.  Values  of  ?4f  up  to  3000  kg/em2  can  be  achieved 
depending  on  the  initial  pressure  ratio  across  the  piston. 

This  limitation  is  imposed  by  the  strength  of  the  barrel  and 
its  end  fittings. 


Leybold  vacuum  diffusion  pump  (type  PD  1000/635)  and  a 
rotary  "backing”  pump  (Type  D12  1?6/l6). 
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FIG.  1 Schematic  of  gas  compression  system. 
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The  pressure  {P4),  temperature  (T4)  and  density  (p4) 
variations  of  the  nitrogen  are  measured  during  the  first 
compression  stroke.  Typically,  the  initial  barrel  pressure 
measured  by  means  of  either  a mercury  manometer  or  a Wallace 
and  Tiefcnan  gauge  is  1 kg/cm2;  the  transit  time  of 
the  piston  from  its  initial  rest  position  to  its  "peak 
pressure"  rest  position  is  about  50  msec.  The  total  measur- 
ement time  of  interest  is  of  the  order  of  1 msec.  During 
this  time,  the  piBton  is  between  10  and  1*0  mm  from  the  closed 
end  of  the  barrel. 

The  facility  and  its  instrumentation  are  shown 

in  Fig. 2. 


2.2  Pressure  instrumentation 


The  pressure  variation  of  the  test  gas  during  the 
compression  is  measured  using  either  a Type  6221  or 
Type  6201  Kistler  piezo-electric  pressure  transducer  and 
a Kistler  568/M5  charge  amplifier  and  recorded  on  a 
Tektronix  Type  502A  Oscilloscope.  The  transducer  is  mounted 
in  the  wall  and  near  the  end  of  the  barrel,  as  depicted 
in  Fig. 3.  The  oscilloscope  trace  is  photographed  irith  a 
Polaroid  camera. 


2 . 3 Temperature  instrumentation 

The  temperature  of  the  test  gas  is  determined 
by  using  the  sodium-line  reversal  technique  (Ref. 8).  A single 
source  effective  double  beam  system,  shown  schematically  in 
Fig.it  is  used.  Two  back  silvered  prisms  are  mounted  in 
the  end  wall  of  the  barrel,  farming  a Bmall  enclosure  into 
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DISPLACEMENT  TRANSDUCER 


FIG.  3 Schematic  of  test  section  of  compression  tube 
illustrating  the  instrumentation  position. 


11 


AEDC-TR-74-47 


FIG.  4 Sodium  line  reversal  optics. 
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which  some  of  the  nitrogen  teat  gas  is  compressed  during 
the  compression  stroke  (see  Fig. 3).  Light  from  a tungsten 
strip-filament  lamp  is  modulated  by  a mechanical  chopper 
and  is  focused  by  a lens  system  and  directed  by  one  of  the 
prisms  into  the  enclosure.  The  light  passes  out  of  the 
enclosure  by  means  of  the  second  prism.  A mirror  and  lens 
system  direct  the  light  to  an  EMI  955SC  photomultiplier. 

Suitable  stops  limit  the  solid  angle  and  the  enclosure 
volume  monitored  by  the  photomultiplier.  The  signal  from  the 
photomultiplier  is  recorded  on  the  Becond  channel  of  the 
Tektronix  502A  oscilloscope. 

The  test  gas  is  seeded  with  sodium  atoms  by  placing 
a few  milligrams  of  sodium  azide,  NaN3  , at  the  entrance 
of  the  enclosure.  Sodium  azide  decomposes  exothermically 
at  600°K  to  form  free  sodium  atoms  and  nitrogen  (Ref. 9). 

Hence  during  the  compression  process,  sodium  atoms  are  formed, 
raised  to  a high  temperature,  and  compressed  into  the  small 
enclosure.  At  temperature  above  1200°K,  the  atoms  emit  light 
at  characteristic  wavelengths  and  this  light  is  monitored 
by  the  photomultiplier. 

The  principle  of  the  sodium  line  reversal  technique 
is  briefly  as  follows.  When  light  from  the  lamp  is  present, 
and  provided  that  its  intensity  is  greater  than  the  intensity 

X 

of  the  thermal  emission  from  the  sodium  atoms  , these  atoms 

will  absorb  some  of  the  light  from  the  lamp.  By  comparing 

the  absorption  and  emission  of  the  sodium  atoms  with  the 

intensity  of  the  light  emitted  by  the  lamp,  the  "electronic" 

temperature  T , of  the  sodium  atoms  is  determined  from  the 
s 

statically  calibrated  brightness  temperature  T_ , of  the 

il 

lamp.  The  relevant  equations  and  assumptions  governing 
the  use  of  the  technique  are  presented  by  Lapworth  (Ref. 10). 

X 

This  is  always  the  case  in  the  present  experiments.  If  the 
contrary  occurs,  the  reversal  technique,  as  used  here,  becomes 
a very  inaccurate  thermometer  (Ref. 11). 
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In  the  present  experiments,  the  spectral  region  of 

interest,  namely,  the  sodium  D-line  doublet,  is  isolated  by 

o 

a filter  vhich  has  a half-band  width  of  UO  A centred  about 

o 

a wavelength  of  5095  A.  If  vj  is  the  photomultiplier  output 

voltage  when  the  lamp  light  is  cut  off  by  the  chopper,  v2  ifl 

the  voltage  when  no  sodium  atoms  are  present  and  only  the 

lamp  light  is  monitored  and  V3  is  the  voltage  when  both 

the  light  from  the  lamp  and  the  sodium  atoms  is  recorded, 

then  T may  be  determined  (Ref. 11 ) from 
s 


1 

T 


+ 


v2  - 


V1 


(1) 


where  k is  the  Boltzmann  constant,  h the  Planck  constant  and 
v the  central  frequency  of  the  sodium  D-lines . vj  is  a measure 
of  the  intensity  of  the  thermal  emission  from  the  sodium 
atoms  alone,  v2  is  a measure  of  the  intensity  of  the  lamp 
light  and  v3  is  a measure  of  the  emission  of  the  sodium  atoms 
and  the  absorption  of  these  atoms  relative  to  the  light  from 
the  lamp.  v2  is  recorded  prior  to  the  compression  stroke, 
v3  and  vj  may  be  measured  almost  simultaneously,  provided 
the  lamp  light  chopping  frequency  is  sufficiently  high. 

Hentall  et  all  (Ref. 12)  show  that  the  electronic  temperature 
of  the  sodium  atoms  is  equal  to  the  vibrational  temperature 
of  the  nitrogen. 

Under  the  present  high  pressure,  high  temperature 
circumstance,  the  vibrational  energy  relaxation  time  of  the 
gas  is  of  the  order  of  0.1  ys.  Therefore,  the  test  gas  is 
undoubtedly  in  equilibrium  during  the  compression  cycle 
and  the  equilibrium  temperature  variation  of  the  gas  iB 
determined.  The  line  reversal  technique  aB  used  here  can  be 
shown  to  measure  the  test  gas  temperature  with  an  uncertainty 

of  less  than  two  percent  (Ref.11)  provided  T_  is  always  greater 

L 
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than  T . 
s 


The  Bensit iveness  of  the  measurement  is  illustrated 
by  the  fact  that  at  2500°K  for  the  sodium  vavelength  a factor 
of  tvo  in  light  intensity  corresponds  to  a change  of  only  - 
l60°  in  temperature  i.e.  around  6 per  cent  (Ref. 8). 


2 . k Extension  of  the  range  of  temperature  measurement 
beyond  2500°K 

Some  research  effort  has  been  applied  to  selecting 
an  accurate  method  of  measuring  the  temperatures  above  the 
2,500°K  limit  imposed  by  the  use  of  a tungsten  ribbon  lamp. 

Tvo  methods  are  chosen  to  be  most  appropriate. 

i } Sodium-1  ine^reyer salute chnigue_UBirig_a_carbon_arc_ source 

Null  and  Lozier  (Ref. 13)  have  reviewed  the  use  of 
a carbon  arc  source  as  a radiation  standard.  Positive  electrode 
crater  temperatures  of  3,0OO°K  can  be  achieved  and  maintained 
with  a reproducibility  of  ±10°K.  The  spectral  radiance  appears 
close  to  that  typical  of  a black  body  of  3,800°K  over  most 
of  the  visible  range,  but  particularly  in  the  vicinity  of 
the  sodium  D-line.  The  spectral  eraissivity  of  the  crater 
approaches  unity  over  this  desired  spectral  range.  This  light 
source  enables  the  existing  sodium  live  reversal  optical 
components  to  be  used  to  extend  the  range  to  about  3>500°K. 

This  Bource  is  a constant  temperature  device,  and  filters 
have  to  be  used  to  decrease  the  effective  temperature  to  lower 
values.  A high  pressure  xenon  arc  lamp  was  also  considered, 
however,  considerable  difficulties  are  anticipated  in  its 
calibration  because  its  spectral  distribution  is  not  close 
to  that  of  a black  body. 
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ii)  Infra-red  technique 

In  this  technique,  the  response  of  an  infra-red 
sensor  (such  as  a fast  response  photo-conductive  indium 
antimonide  detector  cooled  by  liquid  nitrogen)  to  a black- 
body  source  radiation  passing  through  the  test  gas  is  used  to 
measure  the  gas  temperature  (Ref. 3).  An  advantage  of  this 
system  over  the  sodium  line  reversal  technique  is  that  the 
introduction  of  an  appropriate  seed  material,  carbon  mono- 
xide, into  the  test  gas  is  simpler  to  control.  Furthermore, 
the  temperature  range  is  much  greater  (less  than  600°K  to 
greater  than  k,000°K  is  claimed  by  Lapvorth  in  Ref.  Ik). 

An  appropriate  black  body  source  is  described  in  Ref. 15. 

An  accuracy  of  2 % up  to  a temperature  of  k,000°K  is  claimed 
for  this  technique,  but  for  low  pressure  conditions. 

The  first  technique  has  been  selected  for  this  work 
for  two  reasons.  No  experience  has  been  obtained  in  the  use 
of  infra-red  techniques  at  the  Institute.  Furthermore,  none 
of  the  equipment  necessary  for  the  technique  exists  here. 

A suitable  carbon-arc  source  manufactured by  Spindler 
and  Hoyer,  Germany  has  been  selected  and  received.  Experience 
in  the  use  of  the  equipment  is  being  amounted  for  the  future 
extension  of  the  work  close  to  3,500°K. 


2 . 5 Density  measurement 

The  density  variation  of  the  compressed  nitrogen 
is  determined  from  accurate  measurements  of  the  piston's 
position  relative  to  the  barrel  during  the  final  stages  of 
the  compression  stroke. 
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Tvo  complementary  systems  are  used.  A small  lead 
pin,  mounted  in  the  end  wall  of  the  barrel,  is  crushed  by  the 
piston.  The  length  of  the  crushed  pin  gives  the  displacement 
L of  the  piston  from  the  end  wall  at  peak  pressure.  Away 
from  peak  pressure,  the  piston's  position  is  monitored  at 
intervals  of  one  millimeter  by  a displacement  transducer 
mounted  in  the  side  of  the  barrel.  The  pin  and  transducer 
positions  are  shown  in  Fig. 3.  The  latter's  position  prevents 
the  transducer  from  being  subjected  to  pressure  levels  above 
30  kg/cm2.  The  transducer,  a commercially  available  Vibrometer 
TW6- TOO/A  eddy-current  contactless  displacement  transducer, 
monitors  accurately  machined  grooves  on  an  aluminium  sleeve 
attached  to  the  piston.  The  transducer  is  excited  by  a Vibro- 
meter type  100  TRI/A  carrier  amplifier.  The  operation  of  the 
transducer  is  illustrated  in  Fig. 5* 

Knowing  the  piston  position  relative  to  the  trans- 
ducer, the  displacement  (L)  and  volume  (V4)  of  the  compressed 
gas  may  be  determined  from  the  geometry  of  the  barrel  and 
piston.  The  number  of  moles  H,  of  gaB  in  this  volume  is  known 
from  a measurement  of  the  initial  volume  V4^  , pressure  P4^, 
and  temperature  T4^.  The  density  of  the  gas  is  given  by 


P 4 


M 

V4 


v4. 

RT4.  V4 


(2) 


where  R is  the  gas  constant  per  mole.  In  deriving  Eq.2,  it  is 
assumed  that  no  gas  leaks  past  the  piston  during  the  compres- 
sion stroke. 
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2 . 6 Compression  cycle  measurement  control  system 

During  the  final  stages  of  the  piBton  compression 
stroke,  the  pressure,  temperature  and  piston  position  signals 
are  recorded  simultaneously.  The  output  from  the  position 
transducer  triggers  the  sweep  of  the  oscilloscopes.  The 
vertical  amplifier  of  the  Tektronix  type  531A  oscilloscope 
used  to  record  the  position  measurement  operates  in  the 
d-c  mode  and  is  adjusted  so  that  a signal  is  recorded  only 
when  the  grooves  on  the  piston  are  traversing  the  transducer. 
The  time  scale  of  the  trace  is  extended  by  using  the  four 
alternate  trace  capability  of  the  oscilloscope's  type  M plug- 
in unit.  A four  trace  "raster"  display,  which  has  a calibrated 
sweep-back  time  of  30  vs,  is  obtained.  The  oscilloscope 
monitoring  pressure  and  temperature  signals  operates  in  the 
single-sweep  mode  to  prevent  recording  compression  strokes 
other  than  the  first.  Base  lines  on  the  traces  are  added 
just  before  each  test.  To  correlate  the  time  scales  of  the 
two  oscillo8cop«s , a one-kHa  square  wave  of  small  amplitude 
is  superimposed  upon  the  pressure  measurement  and  piston 
position  measurements  signals. 


2.7  Typical  results 

A compression  test  was  carried  out  with  the 
following  initial  conditions  : Pq  “ 30  kg/cm2,  Tq  = 292° K, 

P4.  = 1,285  kg/cm2  and  T_  = 2,155°K.  The  oscillograms  obtained 
from  this  test  are  shown  in  Figs  6 and  7 , On  the  temperature 
and  pressure  oscillogram,  typical  values  of  V| , V2  and  vj  are 
indicated,  together  with  the  photo-multiplier  signal  (v  = 0) 
corresponding  to  the  zero  light  level.  Point  B on  the  position 
measurement  oscillogram  corresponds  to  the  time  when  the  centre 
of  the  first  "tooth"  on  the  piston  arrived  at  the  centre  of 
the  transducer.  At  this  time,  the  distance  L from  the  front 
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Horizontal  sensitivity  = 250  ps  per  square 

Pressure,  vertical  sensitivity  = 525  kg/cn2  per  square 

Temperature,  vertical  sensitivity  = 0,5V  per  square 

FIG,  6 Oscillogram  of  pressure  and  temperature. 


Horizontal  sensitivity  = 5C  ps  per  square 
Vertical  sensitivity  = IV  per  square 

1 - First  trace  3 - Third  trace 

2 - Second  Trace  - Fourth  trace 

FIG.  7 Oscillogram  of  position. 
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face  of  the  piston  to  the  end  vail  of  the  barrel  vas  35-5  mm. 
Each  subsequent  "maximum"  and  "minimum"  corresponds  to  the 
piston  moving  one  millimeter  closer  to  the  barrel  end  vail. 

The  piston  rest  position  at  peak  pressure  is  indicated  by  point 
G on  the  oscillogram.  L is  there  given  by  the  length  (10.59  mm} 
of  the  crushed  lead  pin.  After  point  G,  the  maxima  and  minima 
on  the  trace  give  the  piston's  position  at  intervals  of  one 
millimeter  as  the  piston  moves  avay  from  the  end  vail.  Point  G 
on  the  trace  is  determined  from  the  centre  of  symmetry  of  the 
position  transducer  signal.  Typical  time  marks  from  the  one- 
kHz  square  wave  , superimposed  on  the  pressure  and  position 
measurement  signals,  are  indicated  by  point  tR  on  each  oscillo- 
gram. At  these  points,  an  exact  correlation  in  time,  between  the 
oscillogram  traces,  is  obtained.  Avay  from  these  timeB,  the 
signals  are  correlated  by  using  the  accurately  calibrated  time 
scales  of  the  oscilloscopes. 
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3 . COMPUTER  PROGRAM  AIDS  TO  THE  EXPERIMENTS 

3 . 1 Piston  compression  cycle  •predictions 

Calculations  concerning  the  compression  of  the  gas 
have  been  carried  out  using  a piston  cycle  program.  The  basic 
program  uses  a Uth  order  Runge-Kutta  numerical  method  to  solve 
the  system  of  two  ordinary  differential  equations  resulting 
from  Newton's  law  describing  the  forces  on  the  piston, 
including  friction  during  the  cycle.  The  program  has  been 
modified  in  the  following  several  ways  providing  aids  to 
the  estimation  of  various  sources  of  uncertainty. 

The  gas  in  front  of  and  behind  the  piston  is 
assumed  to  be  in  an  isentropic  state  during  the  cycle.  This 
assumption  is  possible  if  no  shock  waves  are  generated  in  the 
flow.  The  pressure  traces  exhibit  no  sudden  jumps  indicative 
of  shock  waves.  The  accuracy  of  this  assumption  was  examined 
by  modifying  the  program  to  calculate  the  heat  losses  from 
the  gas  sample  during  the  test  assuming  that  heat  conduction 
was  the  only  mechanism  of  heat  transfer.  The  method  of  Knoos 
(Ref.l6)  was  used.  Results  of  this  program  are  discussed  in 
section  5*3. 


The  pressure  of  the  gas  behind  the  piston  was  always 
low  enough  that  an  assumption  of  a perfect  gas  could  be  made. 

It  was  verified  also  that  the  flow  behind  the  piston  could  be 
assumed  to  have  infinite  speed  of  sound  provided  that  the 
piston  speed  never  exceeded  about  on  third  of  local  sound 
speed.  Several  calculations  carried  out  using  the  full  charac- 
teristics solution  behind  the  piston  showed  that  less  than  1 % 
departures  are  incurred  by  this  latter  assumption. 

The  equation  of  state  used  to  calculate  the 
behaviour  of  the  test  gas  in  front  of  the  piston  was  the 
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Enkenhus-Culotta  (Ref.U)  equation  for  nitrogen  and  the  Grabau- 
Brakinsky  (Ref. 17)  tableB  for  air.  Again  an  assumption  of 
infinite  speed  of  sound  was  used.  This  assumption  vas  examined 
by  modifying  the  program  to  calculate  the  behaviour  of  the 
compressed  gas  using  the  full  characteristics  solution,  but 
assuming  a perfect  gaB . Computer  limitations  prevented 
carrying  out  this  calculation  using  a valid  real  gas  equation 
of  state,  but  the  conclusions  of  the  estimations  would  not  be 
expected  to  be  altered  by  using  the  simpler  equation  of  state 
model.  Again  it  vas  shovn  that  1 e b b than  \%  departures  vere 
obtained  assuming  infinite  speed  of  sound  in  the  compressed 
test  gas  for  lov  piston  speeds. 

The  piston  friction  provides  an  important  practical 
consideration  in  evaluating  the  performance  of  the  facility. 

It  is  not  however  an  important  consideration  concerning  the 
measurements  of  the  state  of  the  gaB  Bince  all  three  para- 
meters, p,  p,  T which  define  the  thermodynamic  state  are 
measured  simultaneously.  The  chief  source  of  friction  is 
expected  to  arise  from  the  flexible  piston  seals  used  to 
ensure  no  leakages . A feature  of  the  experiments  to  generate 
equation  of  state  data  is1  that  leakages  provide  a source  of 
experimental  error  since  the  time  dependent  density  variation 
is  calculated  not  only  from  the  measured  piston  motion  but  also 
from  an  assumption  of  constant  mass  of  test  gas  during  the 
cycle.  In  the  piston  cycle  program,  it  was  assumed  that  the 
friction  was  independent  of  piston  velocity,  but  is  a function 
of  both  the  pressure  in  front  and  behind  the  piston.  These 
pressures  act  on  the  flexible  seal.  Various  values  of  friction 
coefficients  were  examined.  The  consideration -of  friction 
causes  a slight  asymmetry  in  the  thermodynamic  measurements 
during  the  compression  and  the  expansion  caused  by  the  rebound 
of  the  piston. 
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The  various  forms  of  the  program  were  used  to  aid  the 
carrying  out  of  the  experiments  as  well  as  to  make  assessments 
of  uncertainties.  Results  from  the  program  are  mentioned  in 
later  sections  of  the  report. 


3 . 2 Data  reduction  programs 

Two  main  data  reduction  programs  have  been  designed. 
One  program  uses  inputs  of  the  measurements  of  the  traces  as 
read  directly  from  the  measuring  table  and  all  the  carefully 
determined  calibration  constants  to  calculate  the  raw  data. 

The  second  program  presents  the  results  to  allow  comparison 
with  appropriate  equations  of  state. 

1 . Data_reduction_of_rav_data 

The  values  of  the  state  parameters  are  generated  by 
introducing  into  the  program  the  co-ordinates  of  the  pressure , 
emission,  absorption  and  position  traces.  The  first  three  are 
introduced  at  equal  time  intervals;  the  latter  is  introduced 
as  the  time,  arbitrarily  referenced,  for  the  distance  trans- 
ducer to  sense  the  equal  interval  distances  defined  by  the 
centre  of  the  teeth  and  the  grooves  of  the  piston. 

The  following  information, the  acquisition  of  which 
is  described  in  a following  section, is  also  provided  : 

- oscilloscope  and  transducer  calibration  curves  and  values, 

- lamp  voltage  and  current,  and  lamp  characteristic  curve, 

- transmission  coefficients  of  the  optical  system, 

- crushed  lead  pin  length, 

- corrections  for  barrel  expansion  and  piston  compression  due 
to  the  high  pressures  developed. 
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The  data  is  then  re-arranged  in  equal  interval 
time  stepB  by  using  a four-point  interpolation  formula. 

A typical  computer  read  out  of  the  program  is  shown  in  Table  1, 


2.  Comgarison_between_experimental_data_and_ theoretical 

state 

A computer  program  was  devised  to  aid  the 
comparison  of  theoretically  determined  equations  of  state  vith 
the  experimental  results.  It  calculates  for  each  measured  point 
(P,  p,  T)  the  following  parameters 

calculated  pressure  from  measured  density  and  measured 

temperature;  p , (P 

calc  meas  r seas 

Also  the  entropies  Si/R{p  , T 

meas  meas 

S3 /R(P , p ) are  calculated. 


T ) i 
meas 


. T ) i “4  Tcalc(P »---)• 


meas 
>■  S2/E(Pmeas- 


meas 

T ) , 
meas 


• P 

meas  meas 


The  theoretical  thermodynamic  information  used  for 
nitrogen  was  the  Enkenhus-Culotta  equation  of  state  (Ref.lt) 
which  agrees  vith  the  AEDC  tables  (Ref. 18)  to  within  a fev 
per  cent  for  all  parameters  over  the  range  considered.  Two 
subroutines  give  pressure  and  dimensionless  entropy  as  a 
function  of  density  and  temperature.  For  air,  the  AEDC  tables 
(Ref. 17)  were  used.  A four  point  interpolation  formula  gives 
pressure  and  dimensionless  entropy  as  a function  of  density 
and  temperature.  In  each  case  a Nevton-Raphson  -searching  proce- 
dure is  used  to  calculate  px,  Tj , S2/R  and  S3/R  when  the 
arguments  are  not  density  and  temperature. 


The  overall  uncertainty  percent  is  also  given,  A typical 
computer  read  out  is  given  in  Table  2. 
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TABLE  I.  Computer  Read-Out  of  Programme  1 (Data  reduction) 


RUN  NUMBER  39GAS#AIR 


TIME  +0R-  RHO  +0R-  LAMBDA  POSITION 
MU SEC  AMGT  MM 


RUM  MUMB.ER  39GAS#AIR 


TIME  +0R-  PRES  +0R-  TEMP  +0R- 
MUSEC  KG 7 CM 2 K 


- 1090,3, 

11. 

4 

22  5.5, 

3. 

1481.2 

13. 

-*10  39,5, 

10. 

4 

240.4, 

3. 

1510. 2 

14. 

-983.7, 

9. 

4 

258.1 

3. 

1543.7. 

14. 

-337,8. 

9. 

i 

2 73.0, 

3. 

1577.6, 

15. 

-88  7.C. 

a. 

4 

291.7, 

1619.3, 

16. 

-836.2, 

8. 

4 

314.0, 

3. 

1648.2 

17. 

-785.3, 

7. 

338.2 

3. 

1665.2 

17. 

-731*. 5, 

7. 

t 

3 58.6, 

3. 

1G72.4, 

17. 

-683.7, 

6. 

4 

385.6, 

3. 

1G8G.0, 

17. 

-G32. 8, 

6. 

4 

413.4, 

3. 

1703.  S, 

18. 

-582.9, 

5. 

l 

443.1 

3. 

1737.4, 

19. 

-531.2, 

5. 

/ 

472. a 

3. 

1735.3. 

19. 

-489.3, 

4. 

4 

509.8, 

3. 

1704.9, 

18. 

-429.5, 

4. 

4 

546.9, 

3 . 

1692. 2 

17. 

-378.7, 

3. 

i 

5 79.4, 

3. 

17C4.Z 

18. 

-327.8. 

3. 

i 

618.3, 

3. 

1731.1 

18. 

-277.9, 

2. 

t 

055.3, 

3. 

1798.3, 

21. 

-226.2, 

2. 

4 

692.4, 

3. 

1815. G, 

22. 

-175.3, 

1. 

4 

721.9, 

3. 

1826.1 

30. 

-124.5, 

1. 

4 

7 5C.6, 

3. 

1836.7, 

14. 

-73.7, 

o. 

t 

775.5, 

3. 

1844. 2 

17. 

-22.8, 

0. 

4 

789.4, 

3. 

1350.4, 

18. 

-1413.3, 

14. 

4 

3 7.5, 

0.1, 

44. a 

44 .22 

-1387.7, 

13. 

4 

38.4, 

0.1, 

45. a 

43.22 

-1353.2, 

13. 

4 

39.3, 

0.2 

46. a 

42.22 

-1317.7, 

13. 

4 

40.2 

0.1, 

47.2 

41. 12 

-1289.4, 

12. 

i 

41.2, 

0.2 

48.2 

40.22 

-1255.2 

12. 

i 

42.2 

0.2 

49.5, 

39.12 

-1223.0, 

12. 

i 

43.3, 

0.2 

50.8, 

38).  22 

-1188.8, 

11. 

4 

44.5, 

0.2 

52. L 

37.22 

-1155.3, 

11. 

4 

4 5.7, 

0.2 

53.5, 

36.22 

-1121.1, 

11. 

4 

47.0, 

0.2 

55.0, 

35.22 

-1088.2  10. 

4 

48.3, 

0.2 

56.  G, 

34.22 

-10  57.6, 

10. 

4 

49.2 

0.2 

58.3, 

33.22 

-1018.2 

10. 

4 

51.3, 

0.2 

60.2 

32.22 

-964.9, 

9. 

4 

53.0, 

0.2 

62. a 

31. 21 

-931.4, 

9. 

* 

54.7, 

0.2 

64.2 

30. U 

-893.9, 

9. 

4 

56.6, 

0.2 

06.2 

29.21 

-859.0, 

8. 

i 

58.6, 

0.2 

68. & 

28.  U 

-815.0, 

8. 

4 

60.8, 

0.2 

71.2 

27. 11 

-778.8, 

7. 

* 

63.2 

0.2 

73.9, 

26.21 

-735.4, 

7. 

4 

65.6, 

0.1, 

76.9, 

25.21 

-692.3, 

6. 

68.4, 

0.2 

80.2 

24. 21 

-645.7, 

6. 

4 

71.3, 

0.2 

83.6, 

25. 11 

-59G.7, 

6. 

i 

74.6, 

0.2 

87.4, 

22. 11 

-550.3, 

5. 

4 

78.2, 

0.2 

91.7, 

21. 11 

-496.7, 

5. 

4 

32.2 

0.2 

96. 4, 

20.11 

-445.8, 

4. 

4 

86.7, 

0.2, 

101.6, 

19. U 

-365.9, 

3. 

4 

91.8, 

0.2, 

107.5, 

18.21 

-298.8, 

3. 

4 

97.5, 

0.3, 

114.2 

17.  21 

-214.0, 

2. 

i 

104.0, 

0.3, 

121.8, 

16. 21 

-98.9, 

0. 

i 

111.5, 

0.3, 

130. a 

15. U 

78.9, 

0. 

i 

119.7, 

0.3, 

140.2 

14.22 
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TABLE  2.  Computer  Read-Out  of  Programme  2 (Comparison  between 
experiment  and  thermodynamic  data) . 

RUM  NUMBER  39  MVR7207 


p 

R 

T 

P(R/ T) 

!UP,T> 

TCR,P) 

S(R,T) 

S(P/T) 

S ( P , R) 

KG/CM2 

AMGT 

K 

KG/ CM2 

AMGT 

K 

225.5 
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4.  CALIBRATION  PROCEDURE 


4 . 1 Pressure  measuring  equipment 

The  pressure  measured  by  a piezo-electric  sensor 
is  calculated  from 


P = 


Ax-K 


where 


p is  the  pressure  (atm") 

K is  the  overall  oscilloscope  sensitivity 

Ax  is  the  oscilloscope  deflection  measured 
photograph  (mm) 

kC  is  the  amplification  factor  (pC/V) 

P 

a is  the  transducer  sensitivity  (pC/atm) 


(volta/mm ) 
on  the 


The  calibration  curves  of  the  pressure  transducers 
used  and  from  which  the  value  of  a is  obtained  were  supplied 
by  the  manufacturer.  The  curves  were  checked  using  a fCistler 
press*  (Kias  Swiss  type  6901  SH  56679).  The  transducer  to  be 
calibrated  and  a Kistler  high  precision  reference  transducer 
are  mounted  on  the  press  through  similar  adaptors.  The  outputs 
of  the  two  transducers  are  monitored  as  the  percentage 
difference  between  them.  The  output  of  the  reference  gauge 
is  monitored  on  a digital  voltmeter.  The  sensitivity  obtained 
by  this  method  as  compared  to  the  maker's  calibration  of  the  two 
VKI  transducers  is  shown  in  Fig.  8. 


This  was  carried  out  at  the  Fabrique  Nationals*  Herstal, 
Belgium 
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FIG.  8 Kistler  pressure  transducer  calibration. 
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The  performance  of  the  charge  amplifiers  used  in  the 

compression  tube  vas  also  checked  out  using  the  press  combined 

vith  tvo  calibrated  VKI  gauges  and  a calibrated  charge  amplifier 

supplied  by  FN.  A fixed  value  of  the  amplification  factor  is 

initially  set  on  the  charge  amplifier.  As  the  pressure  on  the 

transducers  vas  increased,  the  output  voltage  from  the  charge 

amplifier  under  calibration  is  monitored  on  a digital  voltmeter, 

and  the  pressure  is  assessed  from  the  instrumentation  chain 

involving  a calibrated  transducer  and  the  calibrated  charge 

amplifier.  Knowing  the  pressure,  transducer  sensitivity  and  the 

voltage,  the  correct  value  of  k C can  be  obtained  for 

P' 

comparison  vith  the  amplifier  setting.  The  percentage  error 
for  tvo  values  of  k used  in  the  compression  tube  experiments 
is  shown  in  Fig. 9. 

It  is  seen  that  in  both  the  transducer  and  charge 
amplifier  calibrations,  errors  of  the  order  of  would  be 
incurred  if  the  makers  calibrations  were  used.  Corrections  were 
made  in  the  data  reduction  programs  for  these  errors. 


k . 2 Lamp  and  optical  system 

The  Philips  D-8,  6 volt,  16-17  amp.  rating  tungsten 

ribbon  lamp  is  calibrated  against  an  optical  pyrometer  (of 

Evershed  and  Vignoies  Ltd,  London  WU).  For  a black  body 

filament  temperature  Tj,  at  a lamp  current  I and  voltage  V, 

the  effective  brightness  temperature  T at  a frequency  v is 
• ^ vp  P 

given  by  : 


An 

T 


B 
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where  is  the  emissivity  of  tungsten  at  the  frequency 
v and  is  taken  to  he  0.43  (from  tables) 

r 

tj,  is  the  transmissivity  of  the  glasB  envelope  of  the  lamp 

tr  is  the  transmissivity  of  the  prism  and  lenses 

h is  the  Planck  constant  and 

k is  the  Boltzmann  constant 

The  pyrometer  measures  the  brightness  temperature  at  a 

D 

wavelength  of  66 50  A,  whilst  the  experiments  are  carried 

o 

out  at  the  sodium  D-line  wavelength  of  5895  A.  It  is  assumed 

that  there  is  a negligible  error  if.  r , t_  and  e are  taken 

t EFT 

as  constant  for  this  wavelength  change. 


4.2.1  Ejtimat ion  jpf  _t ransm is s_iv it y_  _of_  th_e  J;ungsten_la.iap_gla.s^ 

_e  n_ve_l_o  jr 

E 

When  just  the  lamp  is  viewed  by  the  pyrometer  then  its 
reading  ( ^ 1 given  by 

1 Ik.  k 

It— T * Tl  - ST  ,n  te  ’ iTT  ln  et 

vp  1 F p p 

If  another  similar  lamp,  but  not  illuminated  is  placed 

between  the  first  lamp  and  the  pyrometer,  the  light  emitted 

from  the  tungsten  ribbon  crosses  three  thicknesses  of  glass. 

The  new  reading  of  the  pyrometer  (T  is  given  by  : 

v p 


r 

2 


3k_ 

% 


In 


£n  e 

T 
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The  transmissivity  of  the  envelope  can  he 
found  from  subtracting  these  tvo  equations  to  give  : 


exp 


-hv 

2K 


< 


(T7"r 

p 1 


P 2 


) 


4,2.2  E_stimahio_n  _of  _t r an_s m_i B_s_i vi ty_o_f  _t_h e_  grj. sjns_  _a nd _le nses^  _tr 

This  is  estimated  by  viewing  the  radiating  tungsten 
ribbon  through  the  prism  and  lenses.  The  reading  of  the 

pyrometer  (T  )3  is  given  by  : 

VP 

(T  ) = t7  “ Zn  tE  ■ h^T  ln  et  " hT“  ln  tr 

vp  3 F P y P P 

which  can  be  solved  for  t 

n 

r hv  3 1 

tr  = !iuir  (tt— r - Tr-r1 

up  1 vp  3 J 

Calibrations  carried  out  show  that  both  and 

E R 

vary  with  temperature  as  Bhovn  in  Fig, 10. 


4.2.3  JEstimat ion  _of_  Jthe_e^^ct_i ve_b_right nes s_  _t emp_e_r at ur_e  _at_ 

3 _p a r_t icula r_  fre qjign_£y_ 

The  temperature  of  the  lamp  at  the  centre  of  the  sodium 
D~line  frequency,  v,  is  given  by 
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In  t. 


CT 


T 


R 


The  unknown  T can  "be  eliminated  by  subtraction  from 

F 

the  previous  equation  to  give  : 


1 1 k / 1 1 * . 

— * ~rz \ — - — - — ) in  Em  t 

T ( T„  ) h v v E T R 

v v 3 p 


The  results  of  the  calibration  of  the  effective 
brightness  temperature  of  the  tungsten  ribbon  lamp  used  in  the 
experiments  at  frequency  v as  seen  in  the  test  section  is 
presented  in  Fig, 11. 

The  tungsten  ribbon  lamp  is  povered  by  a stabilized 
6V  16A  power  supply.  The  50  Hz  ripple  of  light  intensity 
due  to  mains  supply  is  less  than  1%  during  the  tests  and  introduce 
an  error  of  0.5%  on  the  measurement  of  gas  temperature. 


4.3  Recording  equipment 


4.3,1  Time_bases 

In  these  experiments  the  most  important  feature  1b 
to  provide  accurate  synchronisation  of  the  traceB  taken  rather 
than  knowledge  of  the  absolute  value  of  the  time  base.  The 
method  of  synchronisation  is  carried  out  by  injecting  on  one 
trace  of  each  oscilloscope  (pressure  and  piston  motion  traces) 
a signal  from  a 1 kHz  square  wave  generator. 


4.3.2  Amplit ude _a ens i t i vi ty 

The  vertical  sensitivity  of  the  oscilloscopes  was 
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FIG.  11  Lamp  calibration. 
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determined  by  injecting  a calibrated  spare  wave  onto  the  beam 
of  the  scope  which  was  photographed  with  the  same  polaroid 
camera  as  used  in  the  experiments.  Both  experimental  and  cali- 
bration traces  were  read  at  the  same  position  of  a trace-reader 
constructed  at  the  Institute  which  itself  was  checked  for  linear- 
ity and  repeatability  by  measuring  many  times  an  accurate  finely 
divided  scale.  The  resolution  of  the  measuring  table  is  about 
5 um.  A cross  wire  on  the  traveller  allowed  easy  selection  of  the 
centre  of  the  oscilloscope  traces  which  were  typically  about 
0*5  millimeters  thick.  Typical  corrections  of  oscilloscope  sen- 
sitivity from  switch  settings  are  of  the  order  of  5% • 

4 . 4 Photomultipliers 

The  EMI  9558c  is  excited  by  a "Fluke"  Type  M412B  high 
voltage  stabilized  d.c.  power  supply  providing  a voltage  range 
from  0 to  2. 1 kV  in  steps  of  0.1  kV.  The  system  is  very  stable 
over  short  periods  of  time,  and  it  is  not  necessary  to  calibrate 
the  system,  since  seconds  before  each  teBt,  a signal  from  the 
calibrated  light  source  is  recorded  on  the  oscilloscope.  The 
measurements  of  temperature  during  the  test  rely  on  this  signal 
as  a reference. 


4 . 5 Position  measurements 

Lewis  et  al  (Ref. 6)  have  described  the  tare  tests 
carried  out  to  ensure  that  the  maximum  output  from  the 
variable  reluctance  transducer  is  obtained  when  the  centre  of  the 
transducer  coincides  with  the  centre  of  a tooth  and  the  minimum 
with  the  centre  of  a groove.  Measurements  using  an  accurate 
travelling  microscope  showed  that  the  grooves  had  been  machined 
such  that  an  accuracy  of  0.2%  over  a span  of  10  grooves  was 
achieved.  It  is  estimated  that  the  crushed  lead  pin  can  be 
measured  to  within  0.1  mm. 
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5.  DISCUSSION  OF  UNCERTAINTIES 


5 . 1 Justification  of  using  a "dynamic"  experiment 
to  measure  thermodynamic  state 

Non-e^uilibrium_behayiour 

During  the  measurement  area  of  interest,  i ..e . the 
last  millesecond  of  the  first  compression,  the  pressure  and 
density  of  the  gas  increase  by  an  order  of  magnitude  and 
the  temperature  is  doubled  (Fig. 12).  As  mentioned  in  2.3  » 
calculations  have  shown  that  over  the  high  density  range 
of  interest,  the  relaxation  time  of  the  internally  excited 
energy  modes  of  nitrogen  and  oxygen  of  the  order  of  0. lysecs 
are  two  to  three  orders  of  magnitude  faster  than  the  rates 
involved  in  the  compression  process,  and  hence  ve  can  conclude 
that  equilibrium  is  maintained  throughout  the  tests  and  the 
errors  arising  due  to  this  assumption  are  negligible. 

Instrumentation  response  provides  another  source  of 
uncertainty.  Provided  that  the  variables  to  be  measured  vary 
in  a smooth  way  as  illustrated  in  Fig. 12  then  instrumentation 
with  a frequency  response  of  1kH2  or  better  would  be  able  to 
follow  the  signals  with  the  predicted  rise  times  with  nfgligable 
attenuation.  All  of  the  elements  in  the  measurement  chain 
have  responses  considerably  better  than  10  kHz.  More  critical 
for.  these  experiments  however  is  the  lag  time  between  signals. 
Although  a lag  of  1 Oysec (typically  equivalent  to  a phase  lag 
of  2°) in  these  experiments  will  cause  only  slight  uncertainties 
near  the  peak  pressure  when  the  pressure  variation  with  time  is 
small,  at  half  the  peak  pressure  the  error  could  be  as  large 
as  3£.  Such  overall  lags  can  be  accrued  from  the  combination 
of  response  times  from  transducers,  carrier  amplifier,  charge 
amplifier,  oscilloscope  amplifiers,  etc.  This  suggests  that 
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only  data  near  the  peak  is  selected  for  study. 


5 . 2 Pressure  variations  in  the  gas  sample 

A theoretical  study  was  initially  made  to  estimate  the 
degree  of  uniformity  in  the  volume  of  the  gas  sample. 

The  compression  cycle  was  calculated  using  the  characteristics 
method  to  solve  the  unsteady  one-dimensional  flow  problem. 

It  was  assumed  that  the  gas  was  perfect.  The  results  show 
that  at  the  lower  pressure  cases  pressure  oscillations 
occur  during  the  compression  but  disappear  during  the  piston 
rebound.  This  is  confirmed  by  experiments,  however,  these 
oscillations  appear  to  be  smaller  than  predicted  (Fig. 13) 


5 • 3 Assessment  of  temperature  variations  in  the  gas 
8 ample 

The  results  of  the  piston  cycle  calculation  outlined 
above  indicate  that  there  are  equivalent  temperature  variations 
through  the  gas  sample.  Since  they  are  connected  to  the  one- 
dimensional,  unsteady  flow  pressure  variations  and  these 
pressure  variations  were  measured  to  be  small,  they  themselves 
can  be  assumed  to  be  small. 

Another  source  of  spatial  temperature  variation  iB 
that  due  to  the  thermal  boundary  layer  on  the  cold  tube  wall. 
The  simplified  piston  cycle  program  described  in  3*1  was  used 
to  determine  the  extent  of  the  thermal  boundary  layer  and  the 
total  heat  loss  assuming  that  the  heat  losses  were  due  to 
conduction.  The  boundary  layer  thickness  was  not  greater  than 
0.02mm  and  the  total  heat  loss  was  not  greater  than  3%  of  the 
total  energy  given  to  the  sample  by  the  piston.  The  losses 
are  not  important  to  the  measurement  since  instantaneous 
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Svveep  time  500  ^.s/D 
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Ap  Scale  02  v /□ 


FIG.  13  Low  pressure  trace. 


AEDC-TR-74-47 


measurements  of  p,  V and  T are  made  during  the  last  part  of 
the  test.  The  sodium  line  reversal  technique  can  be  affected 
by  non-uniformities  in  temperature,  since  the  temperature  is 
averaged  across  the  emitting  and  absorbing  path.  However,  the 
non-uniformity  is  spread  over  only  0.k%  of  the  light  path,  hence 
the  errors  due  to  this  non-uniformity  caused  by  conduction 
heat  losses  vill  not  be  expected  to  be  greater  than  0.4# 


5 • 4 Experimental  study  of  spatial  uniformity  of 
pressure  in  the  test  chamber 

Two  Kistler  type  6221  piezo  electric  pressure 
transducers  vere  placed  in  the  end  wall  of  the  compression  tube 
and  in  the  sidewall  respectively.  The  transducer  in  the  end 
plate  was  flush-mounted  to  the  surface.  The  transducer  on  the 
sidewall,  had  a passage  of  52mm  length  and  2mm  diameter  leading 
to  it.  This  latter  position  is  the  usual  position  of  the 
transducer  during  the  tests  in  which  equation  of  state  data 
is  generated.  Both  outputs  from  the  Kistler  type  568M  charge 
amplifier  units  which  vere  adjusted  to  give  nominally  the 
same  overall  sensitivity  were  monitored  on  an  oscilloscope 
as  well  as  the  difference  between  the  outputs,  but  on  a scale 
fivetimes  more  sensitive.  A series  of  experiments  with  this 
set  up  allowed  an  investigation  to  understand  and  help  decrease 
experimental  uncertainties,  by  making  systematic  changes  and 
observing  the  results.  Frequent  calibrations  with  a Desgra  ges 
and  Huot  deadweight  tester  up  to  500  kg/cm2  vere  made. 

A typical  trace  from  these  experiments  is  also  shown 
in  Fig.  13-  The  output  from  the  trace  monitoring  the  pressure 
difference  which  reached  a maximum  difference  of  ±2 % and  occured 
near  the  peak  was  possibly  due  to  the  following  reasons  : 
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1.  non-linearity  of  the  pressure  transducer  output 

2.  Hon-linearity  of  the  charge  amplifier  output 

3-  Pneumatic  lag  due  to  the  finite  size  of  the  passage 
to  one  of  the  transducers 

4.  Lag  due  to  finite  size  of  diaphragm  thermal  protection 
devices  used  in  some  of  the  tests 

5.  Thermal  response  of  the  transducer 

6.  Frequency  response  of  the  transducer. 

The  first  tvo  considerations  if  dominant  would 
tend  to  cause  deviations  symmetrical  with  the  pressure  peak. 

The  third  consideration  would  tend  to  cause  a deviation 
antisymmetric  to  the  pressure  peak.  The  final  three  considerations 
would  tend  to  cause  asymmetric  responses  but  which  would 
approximately  cancel  out  if  it  is  assumed  the  irregularities 
are  similar  in  each  case. 

The  level  of  difference  and  the  symmetry  of  the 
pressure  difference  trace  shown  in  Fig. 13  tends  to  indicate 
that  non-linearity  of  the  transducers  and  charge  amplifiers 
appears  to  be  the  chief  cause  of  the  difference.  Corrections 
for  non-linearity  as  determined  by  the  calibration  described 
in  section  4.1  were  hence  included  in  the  data  reduction 
of  the  equation  of  state  tests.  Also  ib  provides  evidence  to 
the  calculations  of  section  5.2  which  indicate  that  the  finite 
passage  length  to  the  transducer  causes  negligibly  small 
pneumatic  lag. 

Other  conclusions  obtained  from  similar  tests  are  : 

i)  Thermal  protection  of  the  transducer  diaphragm  prevents 
small  errors  (less  than  1 % ) . There  is  also  no  apparent  lag  or 
decrease  in  frequency  response  by  placing  a thin  asbestos 
pad  on  the  diaphragm  to  protect  it  from  the  hot  gas.  Hence 
it  was  decided  to  mount  asbestos  pad  during  the  experiments. 
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ii)  The  response  of  the  gauge  to  vibrations  is  negligibly 
small 

iii)  The  overall  frequency  response  of  the  measuring  channels 
is  fast  enough  to  cause  only  small  errors  » 

5.5.  Experimental  assessment  of  vail  heat  transfer 
during  a test 

A study  has  been  made  into  the  experimental  uncer- 
tainties of  the  method  of  generation  and  measurement  of  the 
test  slug  associated  vith  heat  transfer  and  thermal  effects. 
Measurements  of  the  heat  transfer  rate  have  been  accomplished 
using  fast  response  thin  film  platinum  resistance  gauges 
mounted  in  the  end  vail.  Many  developmental  difficulties  were 
encountered  in  taking  such  a measurement  in  these  gas  condi- 
tions of  high  pressure  and  high  temperature  and  high  vail 
heat  transfer  rates.  Pressure  sensitivity  of  the  output  of  the 
gauges  (seen  in  early  tests  but  using  a different  design 
of  gauge)  and  vibration  sensitivity  were  measured  to  be 
negligible  by  mounting  the  sensor  in  the  compression  tube, 
alloving  it  to  be  subjected  to  pressure  loads  and  insulating 
it  from  heat  loadB . This  vas  accomplished  in  one  test  by  packing 
its  surface  vith  a 1mm  thick  layer  of  vacuum  greaBe  and  in 
another  test  by  protecting  the  sensor  vith  oil  trapped  in  a 
tube  betveen  the  sensor  and  a close  fitting  but  freely  moving 
piston.  A typical  trace  of  the  direct  output  of  the  gauge  is 
shown  in  Fig.  1 1* . The  heat  transfer  rate  variation  vas 
calculated  from  this  surface  temperature  trace  using  a nume- 
rical solution  similar  to  that  given  in  Ref.  19.  The  value 
of  heat  transfer  rate  thus  measured  is  a factor  of  2 times 
that  predicted  by  the  theory  of  Knoos  (Ref.  1 6)  as  described  in 
Section  3.1  vhich  takes  into  account  heat  losses  only  by 
conduction.  Such  a phenomenon  could  be  caused  by  forced 
convection  of  the  flov  of  gas  caused  by  mixing.  This  mixing 
is  though  to  be  caused  by  a ring-vortex  type  of  motion 
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caused  by  the  oncoming  piston  (Ref. 20),  a flov  configuration 
for  which  it  is  difficult  to  predict  the  heat  transfer  rate. 
This  high  value  of  measured  heat  transfer  rate  provides  an 
important  step  towards  explaining  the  temperature  decays 
that  occur  in  intermittent  facilities  of  the  order  of  8°K  per 
millisecond  at  2000°K. 


5 . 6 Gas  leakages 


Three  circumferential  nylon  flexible  Beals  are 
incorporated  on  the  piston  to  prevent  gas  leakages  (see  Fig. 3). 
One  of  the  seals  is  placed  on  the  rear  of  the  piston  to  prevent 
driver  gas  from  leaking  and  hence  contaminating  the  test  gas 
during  the  piston  acceleration  phase.  The  other  two  seals  are 
placed  in  tandem  close  to  the  front  of  the  piBton  to  prevent 
the  leakage  of  gas  during  the  piston  deceleration  phase. 

The  sealB  are  designed  to  expose  little  of  their  surface 
to  the  test  gas.  0-ring  seals  are  used  everywhere  else  to 
prevent  leaks  from  the  compression  tube.  The  seals  were  tested 
statically  by  loading  the  tube  with  30  kg/cm2  on  one  face 
of  the  piston  and  a vacuum  on  the  other  face.  Only  small 
losses  were  produced  in  times  of  several  hours  on  both 
faces  of  the  piston.  A check  of  the  dynamic  effectiveness 
of  the  sealB  Was  seen  when  a pressure  transducer  was  placed 
in  the  tube  wall  in  such  a position  that  the  piston  passed 
it  during  a test.  A zero  pressure  was  recorded  when  the 
transducer  was  located  between  the  fore  and  aft  seals. 

5 . 7 Contamination  of  the  gas 


5.7.1  Nitrogen 

X 

Dry  nitrogen  with  a purity  of  99*9?  was  used  as 
the  test  gas.  The  impurities  are  mainly  residuals  of  oxygen, 

* supplied  by  Air  Liquide 
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hydrogen  and  water  vapour.  To  minimise  further  contaminants 
being  introduced  by  the  transfer  of  the  gas  from  the  bottle 
to  the  barrel,  the  filling  procedure  was  carried  out  in  the 
following  way. 

The  tube  is  first  cleared  throughout  with  acetone. 
The  barrel  is  then  evacuated  to  the  order  of  ImmHg  pressure 
and  maintained  at  this  condition  for  a few  minutes  to  check  for 
leaks.  The  barrel  is  then  filled  with  nitrogen  to  a pressure 
of  2 kg/cm2,  and  a further  check  for  leaks  carried  out  over 
a few  minutes.  The  tube  is  then  evacuated  once  more  and  the 
barrel  finally  filled  to  the  desired  initial  pressure.  A 
limited  amount  of  contamination  may  occur  during  a test  due 
to  hot  high  pressure  gases  attacking  the  small  exposed  surface 
areas  of  the  O-ring  and  plastic  seals. 

5*7.2  Air 

Ambient  air  was  used  throughout  these  preliminary 
tests  as  the  test  gas.  The  use  of  air,  brought  up  a problem  of 
contaminants  from  the  seals.  Slight  carbon  deposits  were  seen 
on  the  tube  walls  after  high  temperature  tests  , which  are 
assumed  to  come  from  oxidation  of  the  plastic  seals  by  the 
oxygen  in  the  air.  Similar  deposits  were  seen  when  testing 
helium  at  very  high  temperatures  (Ref.21). 

5 • 8 Initial  temperature  of  the  test  gas 

The  initial  temperature  of  the  test  gas  is  taken  as 
the  laboratory  temperature.  This  temperature  is  stable  over 
periods  of  a day  but  varies  between  18°C  and  20°C  during 
the  year.  Corrections  for  this  temperature  to  initial  entropy 


supplied  by  Air  Liquide 
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are  made  in  the  final  data  reduction.  During  the  filling 
of  the  tube,  the  gas  temperature  is  slightly  lowered. 

Since  several  minutes  ellapse  before  the  test  is  carried  out, 
it  is  assumed  that  the  gas  temperature  reaches  that  of  the 
barrel . 


5 . 9 Summary  of  measurement  uncertainties 


The  following  table  summarises  the  results  of  the 
uncertainty  study  of  the  measurements  at  peak  conditions 
presented  in  the  earlier  sub-sections. 


Measurement 

Remarks 

Error 

Pressure 

Transducer  and  charge 
amplifier  errors  assumed 
small  since  directly 
calibrated . 

Oscilloscope  calibration 
and  trace  reading. 

*±1.0? 

Temperature 

25°K  error  in  pyrometer 
reading  at  2,500°K 

50  kHz  ripple  3 % at 
2,500°K 

*±  1 .0J6 
*±0.536 

Position 

Trace  reading  1 % 

Groove  machining  l 0 356 

Lead  pin  measurement  / 
Timing 

*±0.2)6 

*±0 . 356 
m 1 . 0j6 

Synchronisation 

*±0.52 

Overall  error  (as  a sum  of 

above)  * 


Because,  smaller  deflections  avay  from  the  peak 
conditions  are  measured,  and  furthermore  synchronisation  is 
more  difficult,  the  sum  of  the  overall  errors  is  expected  to 
be  in  the  range  of  ±756  to  10/6. 
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6.  RESULTS  AND  DISCUSSION 


6 . 1 Measurements  in  nitrogen 

Table  3 summarises  the  region  in  which  equations  of 
state  tests  were  carried  out  using  nitrogen  as  the  test  gas. 

The  initial  conditions,  and  conditions  defining  the  range  over 
which  pressure,  density  and  temperature  were  measured  are  there 
presented.  The  peak  pressure  at  the  initial  dimensionless  entropy 
for  each  test  are  plotted  on  a Mollier  chart  in  Fig. 15  to 
illustrate  the  overall  range  of  the  tests  and  to  aid  identifi- 
cation of  each  individual  test.  All  14  of  those  tests  in  which 
satisfactorily  readable  oscilloscope  traces  of  all  three  measured 
quantities  are  presented. 


From  the  computer  output  of  the  data  reduction 

programs  (described  in  Section  3-2.2)  are  plotted  p n against 

caic 

p , p , against  p , T , against  T (Fig.l6-lB). 

^meas  calc  & meas*  calc  ° meas  6 

The  deviations  between  the  measured  and  calculated  quantities 
are  plotted  in  Fig. 19  against  an  arbitrary  time.  The  dimension- 
less entropy  calculated  from  each  pair  of  the  three  measured 
quantities  is  plotted  in  Fig. 20.  All  of  the  data  from  all 
fourteen  runs  has  been  reduced  and  plotted  in  the  above 
manner  but  not  included  in  the  report.  The  tabulated  and  plotted 
results  are  held  by  VKI  and  a copy  has  been  forwarded  to  AEDC. 


The  following  discussion  is  based  upon  these  plots  : 

a)  In  some  tests  the  curves  representing  the  compression  and 
expansion  can  be  seen  to  be  not  coincident,  as  one  might 
expect.  This  is  mainly  due  to  errors  in  synchronisation  of 
traces.  The  traces  are  synchronised  with  respect  to  the 
timing  signal,  but  ignores  lags  in  instrumentation,  which  are 
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TABLE  3 Summary  of  nitrogen  tests. 


NO 

P4I 

— 

T41 

P4t 

S4yb 

P(  PEAK) 

P ( MIN  ) 

T C PEAK) 

T(MIN) 

p C PEAK) 

pC  MIN) 

5 

1285 

292.16 

1.159 

22792 

2240 

290 

I860 

1500 

220 

63 

6 

uss 

296.1 S j 

1.147 

2202  8 

73S 

405 

1600 

14*0 

109 

70 

* 

1.206 

296JS 

I.I4J 

22040 

1060 

300 

1705 

1400 

137 

62 

10 

1.000 

2I5JS 

0992 

22090 

1460 

625 

1930 

1600 

154 

75 

12 

1000 

2MJ6 

0.999 

21032 

2520 

550 

2090 

1650 

206 

63 

13 

1.000 

296.15 

Oftflft 

23092 

3020 

300 

2320 

1750 

220 

63 

14 

2000 

297.15 

177  3 

2240ft 

1950 

700 

179S 

1450 

204 

106 

IS 

2.000 

29  ft.  IS 

L767 

22420 

1170 

620 

1570 

1220 

156 

96 

It 

1.2  05 

297.15 

1.139 

22.952 

2920 

530 

2120 

1400 

233 

86 

20 

2000 

297.15 

1.773 

2240ft 

2390 

600 

1830 

1290 

237 

96 

23 

0.700 

296.65 

0.621 

22655 

I960 

300 

2240 

1630 

163 

52 

24 

1.000 

296.95 

as  8 7 

21101 

2530 

450 

2130 

1470 

200 

70 

25 

U)00 

295.15 

0992 

23090 

2720 

340 

2195 

1650 

211 

65 

26 

2000 

29SJS 

1.795 

22.384 

1200 

540 

1550 

1285 

157 

97 
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20  21  22  23  24  29  2S 


INITIAL  PRESSURE  P(ATM)  AT  T=  300‘K 


FIG.  15  Tests  carried  out  with  nitrogen (peat  conditions). 
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CfCALC) 


FIG.  16  P calculated  against  P measured  In 
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TEST  S 


FIG.  18  T calculated  against  T measured  in 
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FIG.  20  Dimensionless  entropy  versus  time  (N£>. 
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difficult  to  predict  to  within  the  10  Msec  accuracy  necessary 
to  keep  the  errors  to  small  values.  A better  method  of 
synchronising  the  signals  is  to  assume  that  the  parameters 
measured,  i.e.  p,  p,  T will  reach  peak  conditions  at  the 
same  time  and  to  re-align  the  signals  accordingly.  This 
assumption  can  be  made  provided  that  the  lags  are  not  so  large 
that  the  signal  is  not  attenuated  by  any  appreciable  amount. 

b)  More  scatter  is  obtained  in  the  measurements  of  temper- 
ature than  in  the  other  two  measurements.  This  measurement 
is  the  least  straight  forward  of  the  techniques  used. 

The  sodium  line  reversal  technique  as  applied  to  the  present 
experiment  gives  accurate  measurement  in  a narrow  region  when 
the  gas  temperature  is  close  to  and  below  the  lamp  temperature 
selected.  Errors  then  get  exponentially  larger  at  lower  gas 
temperatures.  Furthermore  the  chopping  frequency  limits  the 
time  resolution  of  the  curves  of  emission  and  absorption. 

An  interpolation  has  to  be  made  between  the  chopped  traces. 

c)  Borne  of  the  pressure  traces  from  the  low  pressure  teBts 
show  small  oscillatory  deviations  from  the  expected  Bmooth 
variation.  A calculation  of  the  compression  cycle  using  a 
characteristic  solution  for  the  unsteady  one-dimensional  flow 
in  front  of  the  piston  have  shown  that  oscillations  in  both 
pressure  and  temperature  (see  Sections  5-2  and  5.M  exist. 
These  are  not  easily  seen  in  the  temperature  traces  since  the 
light  chopping  is  nearly  in  phase  with  them.  These  pressure 
oscillations  were  smoothed  out  when  taking  the  deflections 
from  the  traces. 

d)  There  is  found  an  asymmetry  about  the  peak  value  in  the 
results  of  temperature  which  is  as  yet  not  fully  explainable. 
This  appears  to  *rise  from  the  distinct  asymmetry  of  the  light 
absorption  measurement  as  illustrated  in  Fig. 5.  The  absorption 
is  seen  to  be  greater  on  the  expansion  part  of  the  cycle  than 
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on  the  compression  part.  The  emission  measurement  however 
is  always  symmetrical  about  the  peak  bb  expected.  The  final 
temperature  is  not  greatly  affected  by  the  absorption  term 
in  Equation  t (Section  2.3),  however,  it  is  felt  that  in 
view  of  the  small  errors  introduced  more  weight  should  be 
given  to  compression  data  rather  than  expansion  data.  The 
following  possible  reasons  are  given  for  this  asymmetry. 

i)  The  optical  system  may  have  been  sensitive  to 
the  vibration  due  to  the  recoil  of  the  facility.  Such 
sensitivity  was  found  to  be  absent  by  examining  the  tempera- 
ture traces  when  the  observation  cavity  was  isolated  from  the 
test  chamber.  There  was  no  discernible  response  to  the  signal 
from  the  chopped  light  beam  shoving  that  the  measuring 
system  was  vibration  insensitive. 

ii)  The  gas  may  become  more  opaque  to  light  the  denser 
it  is,  and  the  more  turbulent  the  state  of  the  gas.  The  latter 
is  expected  to  be  true  during  the  latter  stages  of  the  cycle, 
hence  the  absorption  by  this  process  may  be  more  important 

than  in  the  earlier  stages.  A test  vas  carried  out  after 
paying  particular  attention  to  the  cleaning  of  tube  surfaces 
and  without  introducing  the  sodium  seed  but  no  conclusion  could 
be  made  from  the  results  of  the  test  since  even  the  smallest 
residual  traces  of  sodium  appear  to  remain  emitting  and 
absorbing,  and  a thus  similar  asymmetrical  absorbtion  trace 
was  obtained. 

iii)  The  extra  absorption  in  the  expansion  part  of 
the  cycle  may  be  caused  by  slight  contamination  from  the  sodium 
azide,  ablated  plastic  and  other  sources. 

It  was  not  possible  to  give  any  conclusion  regarding 
this  uncertainty  on  consideration  of  the  results  of  the  test 
described  earlier. 

The  data  obtained  from  the  fourteen  runs  is  reviewed 
by  plotting  all  the  data  in  Figures  21  and  22.  In  the  former 
Figure,  P - has  been  plotted  against  P and  the  latter. 
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T against  T . Figure  21  illustrates  that  there  is  good 

agreement  between  measurement  and  calculation  at  pressures 
below  l400  Kg/cm2,  but  a gradual  trend  for  the  measured 
pressure  to  be  above  that  calculated  at  higher  pressures. 
There  is  a discernable  trend  with  entropy  in  that  the 
deviation  is  slightly  more  marked  for  the  higher  entropy 
cases  than  the  lover  entropy  cases.  Figure  22  illustrates 
that  measured  temperature  gives  good  agreement  with 
calculated  temperature  below  1700°K,  and  a trend  for  the 
measured  temperature  to  lie  below  the  calculated  values  at 
high  temperatures.  In  this  latter  case,  there  appears  to  be 
no  trend  with  different  entropy  eases. 


The  result  outlined  in  this  last  paragraph 
illustrates  that  the  equation  of  state  used  is  predicting 
the  measurements  in  the  region  of  P f 1600  Kg/cm2  and 
T > 1700°K  only  to  within  10  - 15Z.  It  is  noted  that  this 
particular  region  gives  the  largest  deviations  between  the 
AEDC  tables  for  (Ref.l8)  and  the  Enkenhus -Culotta  equation 
of  state  (Ref.U).  However,  no  improvement  in  agreement  is 
obtained  by  using  the  former  theoretical  equation  of  state 
data  presented  in  the  AEDC  tables. 


6. 2 Measurement  in  air 

The  data  from  the  12  tests  carried  out  in  air  have 
been  presented  in  the  same  way  as  for  nitrogen  in  Table  1* 
and  Figures  23  to  30.  In  the  review  graphs  (Figs  29  and  30), 
the  reverse  trend  of  deviation  between  measured  and  calculated 
values  is  found  with  air  than  in  nitrogen  but  the  magnitude 
of  the  deviation  is  considerably  higher  ranging  between  10£ 
and  Uo% . 
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TABLE  4 Summary  of  air  tests. 


NO 

P41 

741 

P 41 

S4I/R 

P(PEAK} 

P(MIN) 

T(PEAK) 

T(MIN  ) 

p (PEAK ) 

p(MlN) 

31 

0952 

29S.1S 

0.85 

23.951 

1498 

344 

20S2 

1746 

200 

37 

39 

0952 

295.15 

095 

23.351 

790 

225 

1950 

14  80 

120 

37 

40 

0952 

296.15 

0.84 

23.961 

727 

390 

1830 

1485 

109 

37 

4! 

0.9  52 

297.15 

0.84 

239.71 

957 

320 

2130 

1660 

135 

37 

42 

0.59 

29795 

0-51 

24.476 

420 

220 

1870 

IS90 

63 

22 

44 

0.59 

299.15 

0.51 

24.489 

680 

240 

1980 

1605 

98 

22 

45 

059 

299.3S 

051 

24.490 

1060 

240 

2110 

1770 

144 

22 

46 

0.55 

300.15 

051 

24.498 

14  30 

390 

2250 

1660 

220 

22 

49 

0S9 

297.75 

0.51 

24.473 

1360 

203 

2260 

1790 

175 

22 

50 

0.9  52 

256.85 

0l84 

23.988 

1295 

417 

22  85 

1770 

203 

37 

51 

0952 

296.95 

0.84 

23.968 

2 505 

433 

2170 

1821 

3 48 

37 

5« 

0952 

295.35 

0.85 

23.953 

1967 

490 

2180 

1710 

266 

37 
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FIG*  23  Tests  carried  out  with  air. 
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FIG.  27  Departures  in  X between  tables  of  air  and 
experiments  versus  time. 
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An  anomaly  occurs  in  temperature  measurement  of 
the  expansion  part  of  the  cycle.  At  times  during  the  expansion 
temperature  is  found  to  he  increasing  whilst  pressure  and 
density  are  decreasing.  A possible  reason  for  this  anomaly 
is  that  the  oxygen  contained  in  the  air  is  undergoing  a 
reaction  with  combustible  components  (such. as  with  the 
plastic  of  the  seals)  thus  changing  the  composition  of  the 
test  gas.  The  results  presented  in  the  Figures  24-30  do  not 
include  data  obtained  during  the  later  stages  of  the  expansion 

Because  of  this  doubt  regarding  the  composition 
of  the  gas,  it  is  difficult  to  say  at  present  whether 
the  deviations  are  due  to  errors  in  the  experiment  or  in 
the  equation  of  state  data  presented  in  the  AEDC  tables  (Ref.  17) 
which  was  used  to  obtain  the  calculated  values  used  in  the  comparisons. 
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7.  CONCLUSIONS 


The  compression  tube  has  been  shown  to  be  capable 
of  generating  samples  of  test  gas  up  to  3000  Kg/cm2  and 
2300°K.  Extensions  up  to  higher  temperatures  are  expected 
to  be  possible.  It  has  been  shown  that  the  selected  instru- 
mentation, piezo-electric  transducers,  sodium  line  reversal 
apparatus  and  eddy-current  distance  transducer  are  suitable 
in  terms  of  response  and  accuracy  to  measure  the  pressure, 
temperature  and  volume  in  the  test  time  of  interest.  Careful 
calibration  has  enabled  overall  uncertainty  I'evels  below 
±52  at  peak  conditions  and  typically  ±72  to  ±102  away  from 
the  peak  to  be  attained.  One  further  uncertainty,  difficult 
to  quantify,  but  considered  to  be  small  for  at  least  the 
nitrogen  tests  is  gas  contamination.  With  further  refi- 
nements, especially  in  synchronisation  of  events  away  from 
the  peak,  the  uncertainty  levels  attained  could  be  reduced 
slightly.  Results  of  tests  carried  out  in  a limited  entropy 
range  in  nitrogen  show  that  agreement  to  within  ±52  with  the 
Enkenhus-Culott a equation  of  state  is  achieved  below  ll»00Kg/cm2 
and  1700  Kg/cm2.  A gradual  deviation  to  values  of  152  is 
found  between  experiment  and  equation  of  state  at  higher 
conditions  in  nitrogen.  The  agreement  between  experiment  and 
equation  of  state  for  air  is  considerably  worse  than  for 
nitrogen  and  the  trends  appear  to  be  reversed.  At  this  stage 
it  has  not  been  verified  whether  these  deviations  are  due 
to  errors  in  published  tables  describing  thermodynamic  state 
or  associated  with  contamination  problems  arising  from  the 
use  of  oxygen  used  in  the  experiment. 
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8.  FUTURE  EXTENSIONS  OF  THE  RESEARCH 

8 .1  Interpretation  of  data. 

A parameter  known  fairly  accurately,  but  as  yet 
unused  in  the  comparison  with  real  gas  equations  of  state  » 
is  the  entropy.  The  compression  is  essentially  isentropic , 
such  that  the  initial  value  of  entropy,  which  can  be  defined 
very  accurately  from  the  initial  running  conditions,  will 
remain  constant  throughout  the  compression  and  expansion.  A 
small  correction  can  be  made  for  the  energy  loss  due  to  heat 
transfer  which  has  been  predicted  and  measured  for  several 
cases.  With  this  redundant  piece  of  information,  more  insight 
can  be  attained  as  to  whether  experimental  deviations  from 
the  equation  of  state  seen  so  far  are  experimental  errors  or 
true  deviations  of  the  equation  of  state  in  defining  nature. 
Further  numerical  routines  must  be  developed  to  examine 
the  results  in  this  light. 

A further  redundant  parameter  whose  use  could  be 
examined  is  the  knowledge  of  energy  during  the  compression 
which  could  come  from  accurate  calculations  and  measurements 
concerning  the  overall  piston  cycle  e.g.  piston  motion, 
friction,  heat  losses,  pressure  variation  behind  the  piston  , 
etc . 


Further  insight  into  the  information  can  be  obtained 
by  fitting  measurements  to  simplified  equations  of  state 
which  give  important  parameters  such  as  the  co-volume  term, 
b,  in  the  Abel-Noble  equation. 


p(“  - b)  - RT 
P 
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and  the  inter-molecular  attractive  force  term,  a,  in  the 
Van  der  Waals  equation 


(p  + ap2 ) (~  - b)  ■ RT 


Also  b,  in  equation  of  Btate  in  virial  form  for 
rigid  sphere  molecules 


- 1 + bp  + 0.625  bp2  + 0.2869  bp 3 + O.155  bp1*  ... 

RT 

Furthermore,  the  parameter  Z»  can  be  calculated 

using 

2 = — 5— 

pRT 

8 . 2 Extension  of  experiments  to  cover  a larger 
temperature  and  entropy  range  in  nitrogen 

Further  experiments  in  nitrogen  will  be  made  to 
extend  the  experiments  of  entropy  to  both  higher  and  lover 
values,  and  to  increasing  the  upper  limit  of  the  temperature 
measurements  to  above  2,500°K  to  as  high  as  possible  belov 
the  limiting  temperature  of  the  carbon  arc  Bource.  This  will 
entail  gaining  experience  in  the  use  of  the  carbon  arc 
light  source  as  a source  for  the  sodium  line  reversal 
technique.  One  result  of  extending  the  range  to  higher 
temperatures  is  that  smaller  volumes  of  gas  are  achieved 
in  the  test  region  and  hence  the  uncertainty  in  p will  be 
increased . 
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8 . 3 Diagnosis  of  the  large  deviations  front  equation 
of  state  or  experiments  in  air 

Experiments  will  be  carried  out  to  attempt  to 
determine  the  cause  of  the  large  deviations  of  the  experimental 
data  taken  with  a test  gas  of  air  compared  to  existing 
equations  of  state.  Manipulations  of  the  data  as  described  in 
Section  8.1  will  be  carried  out  to  aid  finding  the  source 
of  deviation,  whether  experimental  or  true  errors  in 
defining  the  equation  of  state.  A spectroscopic  study  of  the 
test  gas  during  a test  is  recommended  in  order  to  determine 
the  nature  of  any  impurities  that  may  be  present.  If  the 
study  indicates  that  the  deviations  are  due  to  experimental 
errors,  then  appropriate  changes  to  the  experiment  will  be 
made  where  possible  such  as  investigating  the  use  of  soft 
metal  seals  instead  of  plastic  seals.  When  satisfactory 
results  are  obtained  then  similar  extensions  of  the  experiments 
will  be  made  as  for  nitrogen  indicated  in  Section  2. 


8.4  Theoretical  survey  for  an  appropriate  equation 
of  state  for  air 


A literature  survey  should  be  carried  out  to 
determine  an  appropriate  equation  of  state  for  a mixture  of 
gases  such  as  air. 
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